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Abstract 

In Simpson's (1927) classical derivation of the temperature of the Earth in the semi -gray model, the surface temperature diverges 
as the fourth root of the thermal radiation's optical depth. No resolution to this apparent paradox was yet obtained under the strict 
semi-gray approximation. Using this approximation and a simplified approach, we study the saturation of the runaway greenhouse 
effect. 

First we generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls 
on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength 
domain of the absorption. 

Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the 
greenhouse effect saturates and that any Earth-like planet has a maximal temperature which depends on the type of and distance 
to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which 
the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the 
visible, could have at most a surface temperature of 1200-1300K irrespective of the nature of the greenhouse gas. 

We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the 
value of A cut , the wavelength above which the atmosphere becomes optically thick. Unless A cut is small and resides in the optical 
region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. 
This has an interesting observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, 
saturation takes place by radiating through windows in the FIR. 

Keywords: Keywords: Radiative transfer, terrestrial planets, extrasolar planets, Greenhouse 



1. Introduction 

The runaway of the greenhouse effect was first described 
by Simpson ( 1927 T928J. He pointed out to what appeared to 
be a paradox. In his solution, the atmospheric temperature di- 
verges as the fourth root of the optical depth of any greenhouse 
gas, such as water vapor. Without any stabilizing negative feed- 
back, this divergence became to be known as the "Simpson's 
paradox" ( cf.|Goody & Yung|1989| >. 

|Brunt| ( |1929 1 elaborated by combining the effect of heat 
transfer by radiation and turbulence through the assumption that 
both processes can be described with additive diffusion coeffi- 
cients. 

Later, |Sagan| ( |1960| > and |Gold| (1964) noted that if a planet 
is too close to the sun, enough water vapor would evaporate to 
give rise to a positive feedback and heat the surface even more. 
This runaway effect would eventually evaporate the whole body 
of water and increase the atmospheric optical depth in the IR. 



King ( 1963 ) treated the semi -gray atmosphere while implic- 
itly assuming that the frequency separating the incoming short- 
wave (SW) and the outgoing far infrared (FIR) radiation is the 
same as the frequency below which the atmosphere is opaque. 
In addition, King assumed that rsw, the optical depth in the vis- 
ible range, is large as well. King demonstrated that the surface 
temperature diverges and depends through the Chandrasekhar 
//(^-function (Chandrasekhar, 1960) on 77 = t F i R /tsw — the 
ratio between the Far IR optical depth and the shortwave one. 

Later, Shultis & Kaper ( 1969) extended the results of King 
(|1963|l and | Wildt ( |1966fc to include anisotropic scattering. 
Stibbs ( 197l| l extended King's results and obtained the emer- 
gent atmospheric radiation "for any combination of values of 
the parameters". However, a tacit assumption in all the above 
treatments of the problem is that the wavelength where the stel- 
lar and planetary radiations fields are equal, and the wavelength 
at which the strong IR absorption starts, are identical. 

Stibbs defined n — l/r/ as "a measure of the effectiveness 
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with which thermal photons generated by the absorption and 
degrading of the dilute incident visual radiation". In the case 
n <K 1, incident photons can penetrate the atmosphere much 
more easily than the generated thermal photons can escape. 
Consequently, the thermal photons tend to congregate where 
they are generated, except those near the surface, and thereby 
establish a source function in the thermal radiation which in- 
creases with optical depth in the atmosphere to produce the 
greenhouse effect. Interestingly, Stibbs also found the condi- 
tions for limb brightening (when heating at the top of the atmo- 
sphere dominates the structure of the atmosphere), and demon- 
strated that the conditions depend solely on r\. 

|Komabayashi| ( |1967[ ) was worried from the "Possibility of 
exponential boiling-off of the ocean on planets with positive 
feedback leading to endless warming between the vapor pres- 
sure and the greenhouse effect" and attempted to provide a "qual- 
itative analysis using a model of one component-two phase sys- 
tem". Kombayashi was concerned with whether it is possible 
that there is no equilibrium solution to the problem of radiative 
transfer with water vapor leading to a runaway. Komabayashi 
solved the problem by adopting the Eddington approximation 
in the grey approximation, and found a solution in which the 
water vapor did not prevent temperature divergence. 

Ingersoll (1969) discussed the possible runaway on Venus 



and considered only the IR flux from below. He used the stan- 
dard Simpson solution, namely: B*{t, oi ) = F/2jt(2 + (3/2)T for ) , 
where F is the value of the upward radiative flux, T tot is the op- 
tical depth from infinity to the surface, and B(T) = o~T A /n. 
Ingersol discussed the radiative-convective equilibrium model 
coupled to a model in which the IR absorption is due to vapors 
in equilibrium with its liquid or solid phase. This analysis led 
to the conclusion that equilibrium is impossible when the solar 
constant exceeds a critical value as an equilibrium model re- 
quires that the condensed phase evaporates into the atmosphere. 
Ingersoll concluded that a runaway is possible for a planet with 
oceans at the orbit of Venus. 

The works by |Simpson| ( |1927| l, Komabayashi ( j 1967| > and 
Ingersoll ( 1969} gave rise to the so called SKI limit, which is 
the largest outgoing long-wave radiative flux which a planet can 
emit while sustaining the liquid (or solid) phase at the surface. 
When the absorbed stellar radiation surpasses this limt, there 
is no equilibrium and a complete evaporation of the condensed 
phase occurs. 

Rasool & De Bergh ( \ 910) discussed the accumulation of 



CO2 in Venus' atmosphere and its runaway. The authors applied 
the Eddington approximation for Planck mean of k. Thus, the 
absorption was assumed to be constant over all wavelengths, 
and the radiative-convective model for the mean absorption was 
implemented for all wavelengths. 

Kasting ( 1988) discussed questions regarding runaway, moist 



water, but as the solar intensity increased, the limit was reached 
and induced a runaway. 

The increase of stratospheric water vapor with surface tem- 
perature sets an even more stringent limit of 1 .IS©. Thus, the 
liquid water should have evaporated at an earlier time. Mars 
on the other hand, is within the "maximum greenhouse" limit 



(Kasting 1993 1, and remains outside "the first CO2 condensa- 
tion" limit. 



When discussing how the runway takes place, Kasting 



1988) raised the idea that the planetary temperature cannot 



grow indefinitely because there always exists a sufficiently high 
temperature for which the thermal radiations' short wavelength 
tail would eventually be large enough to irradiate the stellar 
short wave radiation reaching the surface. This is one of the 
two primary ideas we elaborate upon here. 



Renno et al. 



( 1994 1 included the hydrological cycle in a 
ID radiative-convective model. The radiation flux was calcu- 
lated using a parametrization scheme and was not solved ex- 



plicitly. Renno ( 1997 1 studied the effect of water vapors and 
showed that there can be two linearly stable solutions to the 
radiative-convective equilibrium when the net forcing is larger 
than a critical value. However, a finite amplitude instability 
can lead to a runaway greenhouse regime when the solar forc- 
ing is larger than a second critical value. The first equilibrium 
corresponds to an optically thin atmosphere while the second 
equilibrium corresponds to an optically thick, highly nonlinear 
atmosphere. Unfortunately, the calculation is highly uncertain 
because clouds were not included in the modeling. 



Nakajima et al. ( 1992 1 performed a study of the runaway 



greenhouse effect with a one-dimensional radiative-convective 
equilibrium model. The publication also supports the point of 
view that "there is no saturation of the greenhouse effect re- 
sulting from an unlimited accumulation of greenhouse gases in 
the atmosphere of a planet." Note that these authors refer to 
the state in which oceans cannot exits in an equilibrium, as the 
runaway greenhouse state. This is not exactly what we discuss 
here. 



Pujol & Fort ( 2002 ) elaborated the SKI analysis by allevi- 



greenhouse atmospheres and the evolution of the Earth and Venus 
and argued that the reason for the Venus's atmosphere runaway 
is the higher (a factor of 1.9 \S @ ) solar insolation, where S® is 
the solar insolation at the mean distance of the Earth from the 
Sun. Using the radiative-convective model, he found that the 
runaway limit is an insolation of 1.4S®. Thus, Venus accord- 
ingly, was in the past below the runaway limit having liquid 



ating the implicit assumption of Simpson ( 1927 1, [Komabayashi 
( |1967| l and Ingersoll] ( |1969| l that the absorbing greenhouse gas 
does so equally at all infrared wavelengths (i.e., gray atmo- 
sphere). The authors found several equilibria and that the anti- 
greenhouse effect can remove the radiation limit. 

The origin of the divergency in the semi-gray model is the 
assumption that radiation can only escape through the optically 
thick IR. If so, an increase in the optical depth requires a higher 
temperature gradient to carry the same heat flux from the sur- 
face. The temperature increase would stop only if heat can 
be carried out by other means. For example, Wea ver & Ra-| 
manathan ( 1995] > extended the Goody and Yung's analysis of 



the Simpson's paradox by mainly considering windows in the 
absorption, that is, they deviated from the original semi-gray 
model. We discuss this possibility and show that it does not 
solve the fundamental paradox. Another important aspect is 
that a high temperature gradient can excite very efficient con- 
vection which will carry the excess heat from the surface and 
limit the temperature gradient to the adiabatic one. We defer 



2 



the discussion of convection to a forthcoming publication. 
Leaving the gray approximation, 



Kasting & Ackerman 



(1986)) used a band model in the diffusion approximation and 
obtained extremely high temperature for Venus, of the order of 
1400K. The diffusion approximation for the radiation field in 



gray atmosphere was extended by Minin (1988) to cases with 
heat sources in the atmosphere. This was carried out by assum- 
ing two B(f) functions, each describing another radiative flux. 
Minin wrote the general equation, though he did not refer to 
King's results. 

Later, |Lorenz et al. | ([1999 ) discussed the climate stabil- 
ity of Titan. Their solution for the radiative transfer included 
the blocking of shortwave radiation and the emission through 
windows in the FIR. They employed the analytical model of 
Samuelson (1983) and the more advanced models developed 
by |McKay et al. | ( 1989| l. In principle, the two treatments as- 



sumed the diffusion approximation for the thermal radiation and 
yielded an expression of the form crT^ m . = a (q + 3t/4) + b, 
where T sur is the surface air temperature, while a, b and q are 
constants which depend on Titan's parameters, t is defined as 
the thermal optical depth of the atmosphere. The above result 
is written so as to expose the dependence on the optical depth, 
showing the resemblance to Simpson's result. The same so- 



lution for the temperature was also applied by |McKay et al. 



( 1999 ), who discussed the anti-greenhouse effect on Titan and 
the early Earth, an effect which arises from a large shortwave 
optical depth. This was also discussed by |Pujol & Fort (2002). 

Mis kolczi] ( |2007[ l repeated the calculation of the radiative 
transfer in the diffusion approximation, essentially rederiving 
the results of both Simpson and Goody & Yung| ( |1989[ l for 
r sul (T fo( ). Miskolczi was bothered by the non-physicality of the 
assumption T tot — > oo, while relying on Milne ( 1922 ) in claim- 
ing that the behavior is not physical. However, his revised so- 
lution still contains the problem that the temperature increased 
indefinitely as the optical depth increases mono tonic ally. 

(2008 ) succeeded to obtain 



More recently, |Rutilya et al. 



analytical results for the semi-gray atmospheres while assum- 
ing the following particular assumptions that: (a) The radiation 
field in the spectral region of low absorption coefficient is pure 
solar, (b) The emission of the atmosphere in the short wave- 
length range can be neglected, (c) The thermal emission of the 
surface is negligible in the short wavelength range. We show 
that these assumptions become invalid as the total absorption in 
the long wavelengths increases. 

In the present paper, we do not discuss the effect of minimal 
insolation leading to a runaway. Instead, we study the satura- 
tion as a function of the total optical depth. To this goal, we 
present a wavelength-dependent analytical and a full numerical 
solution of the semi-gray model as well as simplified analyti- 
cal estimates. We start with several definitions and an exposi- 
tion of the model, and then continue with the numerical results. 
These results demonstrate that there is a marked difference be- 
tween a wavelength dependent model and either the analytical 
diffusion or two streams description within a two or few band 
model. In particular, we show that contrary to Simpson's re- 
sult, the surface temperature does not increase indefinitely as 
the total optical depth increases to infinity, but instead it satu- 



rates even in the classical semi-gray approximation. We explain 
the saturation and derive a simple approximation for the satu- 
ration temperature as a function of the model parameters. The 
advantage of the semi-gray model is in its simplicity. As the 
model ignores the details of the absorbers, the obtained results 
are quite general and do not depend on the properties of any 
particular molecule. 

We begin in ^2]by developing the simplest semi-gray model 
which can describe the saturation of the greenhouse effect. In 
§j3] we discuss the conditions in which the Simpson solution is 
obtained. We then show in ^4] how the solution is modified 
when a window in the FIR is introduced. We then describe the 
numerical results to the full semi-gray model in ^6] The anti- 



greenhouse effect is discussed in § 6.2 and the implications of 
the results are discussed in ^7] 

2. The Semi-Gray model 

To understand how saturation arises in the semi-gray model, 
we construct the simplest model which encapsulates the salient 
features. Namely, we assume an optically thin SW band and 
an optically thick FIR band, as described in fig. [T] A cul is the 
wavelength separating between the shorter and longer wave- 
length ranges of the absorption coefficient. In all treatments 
hitherto it was tacitly assumed that A cut = A ra d, where A ra d is 
the wavelength at which 



fh(T t ,A rad ) = I p {T p ,A rad ) = g 



(1) 



in the atmosphere of the planet. Here, / is the attenuation fac- 
tor of the stellar radiation at the planet's orbit, T t is the effective 
temperature of the star which radiates as a black-body, and T p 
is the surface temperature of the planet which also radiates like 
a black-body. /» is the stellar specific intensity at the stellar sur- 
face and I p is the specific intensity emitted by the surfac^] In 
the case of the present Earth with T p = and atmospheric 

composition, g becomes negligible at A — A ra d = 5x 10 4 A, but 
this is a particular case which quickly becomes invalid as T p , 
or T„ or /, change. The approximate equality in the case of 
the Earth, between the wavelength separating the two absorp- 
tion domains of the microscopic physics and the wavelength 
separating the two domains of stellar and planetary dominant 
radiation is accidental, and need not be the general case. Here, 
we alleviate this particular tacitly assumed equality and discuss 
the general case in which A cut + A ra d- 

To solve analytically the radiative transfer through the at- 
mosphere, we shall assume for simplicity that the atmosphere 
includes only scattering (i.e., "infrared clouds"), and that it can 
be described with the two stream approximation. It also implies 
that the temperature structure of the atmosphere is unimportant. 
This will allow for an analytic solution of the radiative transfer, 
since without thermal redistribution, the fraction going through 
different bands remains the same at different heights. The oppo- 
site limit, of pure absorption, would give similar results but not 



'We denote by /'s either upwards or downwards flux, and by F the net flux 
downwards. If a wavelength is given, it is per unit wavelength, otherwise, it is 
integrated over a whole band. 
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Figure 1: A heuristic description of the semi-gray model. The SW band is 
assumed to be optically thin. The FIR is dominated by an optically thick atmo- 
sphere. Note that we do not assume that all the solar flux is in the SW band nor 
that the thermal radiation of the planet is confined to the FIR. 



identical, because of the frequency redistribution when emitting 
at different temperatures. The full numerical solution which fol- 
lows in ^6] solves the more complicated pure absorption prob- 
lem (the algorithm of which is described in Appendix B with 
and without scattering). 

Let t be the optical depth in the FIR, measured from the top 
of the atmosphere downward. If we write 7+(r) as the thermal 
FIR flux towards larger optical depths (downwards) and 7_(t) 
as the flux towards smaller optical depths, then the solutions for 
7 ± (t) under the above approximations are (cf. Appendix B): 

Fhr(t) = 7_(t)-/+(t) = const. (2) 
27 FIR (r) = 7_(t) + 7 + (t) = [7_(t) - 7 + (r)] r + const. 

or by comparing the conditions at the top (t = 0) to the bottom 
(t = T tot ), we obtain 



7-(T tot ) - 7 + (r tot ) 
7- (r to t) +7+ (T tot ) 



7_(0) - 7 + (0), 
[7_(0)-7 + (0)]T tot 
+ [7_(0)+7 + (0)]. 



(3) 



The boundary conditions we have are: 
7+(0) = 7* ;FIR and 7_(r tot ) = I p ,fir, 



(4) 



where 7* fir is the insolation for A > A cut and 7 ; ,fir is the emis- 
sion from the planetary surface at A > A cul . It is generally as- 
sumed that 7* FIR = and that I p FIR is the total thermal emis- 
sion, an approximation which holds for the Earth at T = 288K. 
However, if A cul decreases, these two assumptions are no longer 
justified. 

The last equation is obtained from the assumption of ther- 
mal equilibrium, i.e., the total absorption equals the total emis- 
sion: 



(1 - a)7*, sw +Ia,i = (1 - a)7 p ,sw + I P ,i 



(5) 



where h,sw is me insolation for A < A cut , Ia.i is the atmospheric 
emission towards the surface (at A > A cut ), I p ,sw me planetary 
emission at A < A cu , and a is the albedo at the short wavelengths. 
Our main point is that I Pt sw must be included in the energy 
balance at relatively high surface temperatures. 

Using the two sets of eqs. [3]and|4] the thermal equilibrium 
becomes: 



(7*,; 



SW ' 



'p,SWl 



2(7 a fi R - 7^fir) 

2 + Ttnt 



(6) 



The behavior for T to t —* is as expected, giving direct 
equilibrium between the incoming stellar flux and the outgo- 
ing thermal radiation. However, the "Simpson limit" of large 
optical depths, together with 7*fir = and 7 p sw = 0, gives 
that (1 - a)7+ sw ~ 7 /)> fi R /(t/2). This is an artifact of the two 
stream approximation. In the full radiative transfer solution for 
this limit, one should obtain (1 - 0)7* s W m 7 p F ir/(3t/4). Thus, 
we should modify eq. [6] to be 



(1 - a) [7*,sw - Ip,sw(Tp)] = 



[7p,FiR(7'p) - 7* fir] 
1 + 3T tot /4 



(7) 



This equation is essentially an equation for T p . 

We assume that both the stellar and planetary emissions can 
be adequately described as black bodies. The stellar radiation 
reaching the surface of the planet in the shortwave is given by 



7*,sw - 



(8) 



Here d is the distance of the planet from the star. The factor 4 
arises from the ratio between the total surface area of the planet 
to its cross-section to the stellar radiation assuming fast rota- 
tion. For convenience we work in frequency space, and with 

^cut — c i ^cut ■ 

The long wavelength stellar energy flux at the top of the 
atmosphere is given by 



7*, fir - 



(9) 



The emission of the planet at short wavelengths is given by: 
,v)dv, (10) 



nB(T p ,i 

and at long wavelengths by 

7 p ,fir = f nB(T p , v)dv. 
Jo 



(ID 



If we define / = {R+ld) 2 /A and plug the fluxes into the sur- 
face equilibrium given by eq. [7] we obtain an equation for the 
equilibrium temperature T p as a function of the stellar emission, 
the distance to it, the optical depth in the IR, and of v cut : 



(1-a) [fBiT^v) - B(T p ,v)]dv 



(12) 
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3. The "Simpson" Solution 

To understand the Simpson paradox, we begin by imposing 
the usual approximations in eq. 12 The "Simpson solution" 



is obtained if the planetary SW emission is neglected on the 
l.h.s., while the stellar component in the FIR is neglected on 



the r.h.s. of the governing equation 12 Furthermore, one also 
assumes that v cut is well below the incoming shortwave emis- 
sion but well above the outgoing planetary radiation. Under the 
latter approximation, the integrals give the Stephan-Boltzmann 
emission, namely: C B(T, v)v = crT 4 /n and eq. 



12 



reduces to 



+ 4 ) 4 I d ) 



1/4 



(13) 



This result demonstrates how the planetary temperature di- 
verges as oc Tfof. The main reason for the divergence is the fact 
that irrespective of how high the planetary temperature is, there 
is no route for the planetary radiation to escape except through 
the optically thick atmosphere. 

All other radiative transfer solutions found in the literature 
to overcome the greenhouse runaway involve a window in the 
IR due to water vapors (or any other species). The exception is 
Kasting ( 1988) who suggested that leakage through the SW can 
in principle serve as such a route. 

Besides radiative transfer, another physical mechanism that 
can advect the energy from the surface towards higher up in the 
atmosphere is that of convection. If it arises, it will necessarily 
reduce the equilibrium temperature, and it therefore implies that 
any solution we later find is necessarily an upper limit. Some 
authors concluded that convection may become important in 
a runaway planet, especially if it is dominated by water (e.g., 



|Kasting 1988|l. Oth ers, on the other hand, found that it is not 
(e.g., Matsui & Abe 1986 1. In the only example we know, that 
of Venus, most of the atmosphere is not convective. 

In this respect, is also worth mentioning that that unlike 
convection, the radiative solution is insensitive to the physical 
length scale. This implies that a given radiative solution may 
be adequate to describe one planet, but convection may appear 
on another with a similar atmosphere and energy budget. 

4. Adding a window in the IR 

Although Simpson realized the existence of a paradox, he 



argued that the water vapor window should eliminate it ( |Simp- 
1928) > but without evaluating the impact of this assump- 



son 



tion on his paradox. Wea ver & Ramanathan| ( |1995| l consid- 
ered the effect of such a window in mitigating the runaway 
and even harnessing it. They applied the semi-gray model, 
as assumed by Simpson, and added a transparent window to 
the old Schwarzschild problem. They allowed the scale height 
of the trace greenhouse gas to be different from the one of 
the main constituents and accounted for strong non-overlapping 
lines within the context of semi-gray optical depth definition. 
Weaver and Ramanathan assumed that the short wavelengths 



range is completely transparent and solved for the IR flux only. 
Adopting the two stream approximation they wrote that: 

= D(I + - nB) and — = -D(F - nB), (14) 
dr dr 

where D is a constant that arises from the integration over a 
hemisphere, and r is the constant optical depth in the IR. Usu- 
ally D = 3/2. 

For radiative equilibrium, the authors imposed the condi- 
tion: dF/dT - 0, where F = I + — /_ = —Fq. Here Fg is in- 
coming solar radiation reaching the surface, which has to be 
radiated upwards in the FIR, thereby giving a negative net flux 
F. For global mean condition for Fq is given by 



F = -(l-a), 



(15) 



where S is the solar insolation. The boundary conditions are: 
/ + (r = 0) = and 7_(t) = nB„ = 7 + (r) + F . (16) 



The solution for the temperature is: 
(l+Dr) 



and 



4 (2 + Dr tot ) 
°~ T » = = ■ 



(17) 



Note that there is a discontinuity between the surface and the 
atmospheric temperature. 

What happens when the IR domain contains a window? 
Weaver and Ramanathan defined a parameter B as follows: 



I h 2 Bi(T)dT 
f Q °° B A (T)dT 



(18) 



where (Ai , A2) is the wavelength interval of the window. 
The solution for the radiation field is then: 

a-T (t) = F - 



'(2 + BDt, oI ) 



(19) 



The surface temperature again has a discontinuity, and it is 
given by: 



t4/ s r (2 + Dt, 0I ) 

crl Jt,„,) — tn . 

pK \2 + BDt w ,) 



(20) 



In the limit of B — > 0, we recover Simpson's solution with its 
paradox, namely, T p — > 00 with t* oi — > 00. However, for a finite 
B, there is no divergence in the limit t w , — > 00. Instead, the 
surface temperature becomes independent of t and tends to: 



lim crTt = —. 

r IO ,^oc p B 



(21) 



According to the authors, B is independent of temperature, 



and the result should be compared with eq. 17 However, this is 
not the case. 

Consider three typical windows in the IR. The first window 
is at wavelength (2x 10 4 A, 3x 10 4 A), the second at (5x 10 4 A, 8x 
10 4 A) and the third at (10 5 A, 1.5 x 10 s A). Their corresponding 
B is depicted in fig. [2] 
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Figure 2: The behavior of p(T) as a function of temperature, for three typical 
windows. The first window is (2x 10 4 A, 3x 10 4 A), the second is at (5x 10 4 A, 8x 
10 4 A) and the third is at (1 X 10 5 A, 1.5 X 10 5 A). 



We find that [3 is not temperature independent. Instead, 
f3(T) first increases with temperature (and lowers the green- 
house temperature) and then decreases towards zero (and raises 
the greenhouse temperature). Hence, we conclude that a win- 
dow cannot resolve the Simpson paradox by preventing the tem- 
perature divergence. 

Pujol & North| ( |2003 1 introduced a semi-gray model into 



radiative-convective ID model with a stratosphere in radiative 
equilibrium and a troposphere fully saturated with water vapors. 
They applied the formalism of Weaver and Ramanathan and 
defined a similar /3. Their eq. (12) is essentially the result of 
Weaver and Ramanathan including the particular assumption 
on the absorption. The authors argued correctly that ft depends 
on temperature. 

The first conclusion reached by these authors was that: "The 
single (absolute) SKI limit found in gray atmospheres is not ob- 
tained in non-gray atmospheres with fully transparent infrared 
regions". In other words "the long-wave radiation emitted by 
any non-gray atmosphere with fully transparent infrared regions 
is not bounded". This is a consequence of not noticing that 
P — > as T p insceases. 

5. An overall view for a scattering atmosphere 

We now turn to solve the semi -gray model (eq. [12^ , in par- 
ticular, while considering that the planet's thermal radiation can 
escape through the SW range, and some of the stellar SW ra- 
diation is absorbed at the top of the atmosphere. Since there 
is no analytic solution to eq. 12 we solve it numerically. The 



results shown here are those obtained for a Main Sequence star 
radiating like a black body at 5800K and a fast rotating planet 
(namely, one with a uniform temperature distribution) placed at 
the orbit of Venus, having a similar albedo. 

Fig.[3]depicts the equilibrium surface temperature as a func- 
tion of A cm for different optical depths. Several points are ev- 
ident. First, for each optical depth, there is an "optimal" A cut 
for which the greenhouse effect is the largest. Second, the equi- 
librium temperature is a monotonically increasing function of 
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Figure 3: The equilibrium surface temperature (eq. 1 12) , as a function of A CM 
for several optical depths, tfir = 0. 1 , 1 0, 1 00, 1 3 from bottom to top. The con- 
ditions are those for a planet which has an albedo of 0.75, located 1.08x 10 13 cm 
from a sun like star. The straight dashed lines denote the values expected from 
the Simpson solution for each optical depth. The latter solution breaks down 
for large optical depths. 

the optical depth. Third, for a given A cut , there is saturation 
of the greenhouse effect. Fourth, a larger A cut > 10 4 A im- 
plies a smaller T tot for which saturation is reached. Finally, for 
A ~ 10 4 A, the largest T tot is needed to reach saturation. 

Fig. [4] presents the equilibrium temperature in the A cul - 
log(rpiR) plane, in which three regions can be identified. For 
low tfir and intermediate values of A cut , the Simpson's solu- 
tion is adequate. In this region, most of the flux leaks out of the 
planet through the moderately thick FIR. Because of this, the 
actual value of A cut is irrelevant. 

If the optical depth becomes sufficiently large, the temper- 
ature increases to such high values that a significant part of the 
flux can leak through the SW range, below A cut . For larger 
/Icut's, in cases that the stellar IR radiation dominates that of the 
planet at even longer wavelengths, it is easier for the radiation to 
leak through the SW (which now contains longer wavelengths), 
and therefore saturation takes place at lower optical depths. In 
the A cm - tfir plane, this takes place in the upper-right region. 
Once leakage through the SW becomes important, the value of 
t F ir is unimportant (provided that it is large enough) since the 
small amount of radiation passing through the optically thick 
FIR is irrelevant 

The temperature decrease on the l.h.s. of the plane with the 
decrease of A cul is a consequence of another effect. If /l cut is 
sufficiently small, the atmosphere blocks some of the incoming 
SW radiation from reaching the surface. Instead, some of the 
stellar radiation will be scattered in the atmosphere and radi- 
ated back to space (through the /_ component in the two stream 
approximation, or more generally through 7^, as in fig.[T|. 

Formally, for very small /t cut , the saturation temperature can 
be quite high, though never more than T t . However, it requires 
extremely large optical depths and it also assumes that there are 
no windows in the FIR. More realistically, at very large wave- 
lengths beyond the rotational molecular lines, that is, beyond 
A > 10 6 A, the atmosphere should be optically thin. To see this 
effect, we plot in fig. [5] the equilibrium temperature obtained 
once we allow the thermal radiation to escape beyond 10 6 A. 
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Figure 4: The equilibrium temperature (eq. |12) as a function of i cut and tfir, 
for Venus like conditions (given in fig. [5}. Black lines denote the planetary tem- 
perature while the gray lines, the ratio between the solution temperature and the 
Simpson temperature. The latter is an adequate solution only for moderate op- 
tical depths and moderate cutoff frequencies. If the optical depth is too large 
(upper right regions), the Simpson solution breaks down because radiation es- 
capes through the optical. If the optical is blocked (left region), less stellar 
radiation can reach the surface, and the Simpson solution breaks down as well. 



Two regions remain the same. The equilibrium temperature in 
the region where most of the radiation escapes through the opti- 
cally thick FIR is determined following the Simpson's solution. 
It is denoted here as region A. In the region denoted here as "re- 
gion B", saturation is achieved through leakage in the SW. In 
the third region, "region C", saturation is achieved by allowing 
radiation to escape through very long wavelengths. This region 
is absent in fig. [4] where there is no long wavelength window. 

One interesting aspect about the latter region is the fact that 
SW radiation is blocked, but radiation can escape through very 
long wavelengths. This gives rise to an anti-greenhouse effect 
(e.g., as is probably present on Titan, McK ay"ef al. |1999] l. 

In fig. [6] we see the saturation surface temperature for a 
planet having a vanishing albedo and located at a distance of 
1.5 x 10 13 cm from a main sequence star, having a surface tem- 
perature of T„ = 5800K. As A cul — > 0, the surface temperature 
tends to the black body temperature of the star. It does not 
diverge irrespective of the dilution factor, as is required by ther- 
modynamics. 



Although eq. 12 does not have an analytic solution, it can be 
solved under the approximation that A cu , resides in the Rayleigh- 
Jeans long wavelenth tail of the stellar radiation, and in the 
short wavelength cutoff of the thermal emission. This is car- 
ried out in Appendix A, where we also show that it typically 
under-predicts the accurate solution by only 20K to 40K. 

6. Numerical results for an absorbing atmosphere 

In addition to the analytic description of the semi-gray prob- 
lem, we can also solve the full problem of radiative transfer, 
while alleviating the assumption of a purely scattering atmo- 
sphere. 
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Figure 5: Similar to fig.[4]with the exception that the atmosphere is assumed to 
be optically thin for A > 10 6 A. The gray regions here define the characteristics 
of the solution. In region (A), more than half of the radiation escapes through 
the optically thick FIR. In region (B), more than half of the radiation escapes 
below /i cut , i.e., in the optical and NIR. In region (C), more than half of the 
radiation escapes through the optically thin window beyond 10 6 A. Note that in 
this region, the system has an anti-greenhouse behavior. This is because SW 
radiation cannot easily penetrate the system, but sufficient radiation can escape 
in the FIR "window". 



The details of the numerical solution for the radiation field 
are given in appendix B. While the optical depth previously re- 
ferred to scattering, from now on it refers to absorption. Also, 
unless otherwise mentioned, we will consider here a Venus like 
planet, i.e., a planet revolving a solar like star at a distance of 
1.08 x 10 13 cm, but with an a = albedo. 

We begin by discussing the tfir — > oo limit, showing again 
that there is a maximal possible temperature for a given A cul . 
We then continue to the general case and study the different 
limits where different physical processes govern the planetary 
temperature. Afterwards, we explore more specific points, the 
effects of a FIR window and the anti-greenhouse effect. 

The result of the numerical solution is shown in fig. [7] where 
the surface temperature is plotted as a function of the total op- 
tical depth in the FIR range, for several total optical depths in 
the visible range. The appearance of saturation for high optical 
depths is clear. As t$w increases, the saturation temperature 
decreases for large Tfir. Moreover, the surface temperatures 
for small tfir's decrease to below the equilibrium temperature 
leading to an anti-greenhouse effect. 

The effect of A cu , is shown in fig. [8] As A cut decreases, the 
saturation temperature increases and reaches a maximum. This 
increase takes place because more of the thermal emission is 
blocked by a high FIR optical depth atmosphere. The subse- 
quent decrease for smaller A cu , takes place because more of the 
stellar insolation is in the optically thick wavelength band. This 
radiation is absorbed at a high altitude and re-emitted to space 
without reaching the surface. Thus, less radiation in the visible 
reaches the surface. 
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Figure 6: The maximum possible greenhouse temperature (eq. |12) , with 
1"FIR — » °° as a function of A cut , for an a = planet at a distance of 1 .5 X 10 13 cm 
from a main sequence star with a surface temperature of 5800K radiating like 
a black-body. In principle, the planetary temperature can asymptotically reach 
the stellar temperature, but this requires A CM to be unrealistically small, and 
tfir to be unrealistically large. 
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Figure 7: The surface temperature as a function of the total optical depth in 
the FIR, calculated for several total optical depths in the visible using the full 
numerical calculation. The model planet is one with a = 0, at the distance of 
Venus from a sun like star, with A cut = 30, 000A. An anti-greenhouse effect is 
obtained for rsw = 10 and tfir < 6. 



6.1. Saturation of the FIR emission 

Let us return to eq. [7] and consider the limit t,„, — > oo. In 
this limit I*,sw = Ip,sw(Tp) irrespective of the albedo. 

Furthermore, by considering that I* t sw - I* - /*,fir (where 
1+ is the bolometric stellar flux) in the energy balance equation 
for the top of the atmosphere, one finds that + I*.fir - 

Ip.S w 



+1a,-\, where Ia,i is the upwards emission by the top of the 



atmosphere, in the FIR range. 
One therefore has in the r tot 



oo limit that 



p,SW 



I 



FIR- 



(22) 



This implies that the emission in the FIR from a runaway green- 
house planet is simply the "reflected" stellar component in the 
FIR, which is less than the value expected for a small r FIR . In 
the latter case, the FIR includes also the "reprocessed" stellar 
SW radiation. This may be somewhat counter intuitive, that 
a planet which has a strong greenhouse actually has a smaller 
emission in the FIR. 
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Figure 8: The variation of the saturation temperature with A cut , as calculated 
with the full numerical calculation for the case described in fig. [7] with rsw = 
0.1. The difference between the t = 500 and t = 1000 cases is too small to be 
discerned, thus demonstrating the saturation. 
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Figure 9: The fraction of the flux coming out in the FIR as a function of the 
total optical depth in the short wavelength range, A < A cut , as obtained with the 
full numerical calculation. 



Let us now study numerically the FIR behavior of a planet 
with a strong greenhouse effect. Fig. [9] depicts the fraction of 
the radiation emitted to space from the top of the atmosphere 
in the FIR, as a function of the total optical depth in the FIR 
range, for several cases of tsw- We observe that as t fir in- 
creases, which also increases the surface temperature, the frac- 
tion of the radiation emitted above the cut-off decreases and 
reaches an asymptotic value. This is consistent with the above 
argumentation. 

Next, let us consider the dependence of /fir/7* as a func- 
tion of T] = kfir/^sw, as shown in fig. 



10 We observe that for 



T] » 1 the emitted radiation obtains an asymptotic value which 
depends on rsw- Moreover, as rj decreases, /fir//* reaches a 
maximum as rj decreases, for r\ ~ 1 . 



We show in fig. 1 1 how the FIR emission varies with the 
surface temperature. The observed behavior is far from triv- 
ial, even the sign of the dependence of /fir with respect to the 
surface temperature changes ! 

6.2. The anti- greenhouse effect 



Let us return to fig. 11 As rsw increases to about unity 



and higher, the surface temperature decreases to below T eq (ob- 
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ent heating. We find that without exception, 
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Figure 10: The same as figure[9]but for A cut = 5 X 10 4 A and as a function of 
V = kftr/ksw- /fir = I(A > -^cut) is the FIR radiation emitted from the star. 
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Figure 1 1 : The same as figure|9] but as a function of the surface temperature. 

tained for vanishing Tsw an d t fir, that is, assuming no atmo- 
sphere, which for the depicted case is about 325K). The reason 
is as follows. When tsw increases to about unity, most of the 
incident radiation is absorbed at low optical depths (high alti- 
tudes) and as a consequence the temperature saturates, or even 
decreases inwards. As tsw increases, less radiation reaches the 
surface and the surface is therefore required to emit in equi- 
librium, less radiation. However, if tfir is smaller than Tsw, 
the smaller required emission from the surface translates into a 
reduced temperature. 

In fact, an optical depth of tsw ~ 10~ 2 is sufficient to 
heat the top of the atmosphere so as to invert the temperature 
gradient. This is the anti-greenhouse effect and it appears for 
tsw ^ 10~ 2 and small tfir's. But as tfir increases to about 
Tsw ar, d more, the classical greenhouse effect (hot surface and 
negative temperature gradient with height) re-emerges. A sim- 



ilar solution was found by McKay et al. ( 1999| l, who used a 
radiation diffusion model. 

It is obvious that there exist for every A cut and tsw a value 
°f Tfir for which the insolation results in an almost isother- 
mal temperature gradient in the atmosphere despite the appar- 



( dT m \ 
\dT SW J 



<0 



(23) 



Tp[R=const 



namely, an increase of the absorption in the short range de- 
creases the surface temperature simply because stellar energy 
is absorbed higher in the atmosphere and can more easily es- 
cape from the planet. The magnitude of the derivative depends 
on tfir and increases (in absolute value) with tfir. 

7. Application to Venus 

In relating the results of the idealized model to a real sys- 
tem, we have to be aware that the relevant choice of A cM and 
tfir depends on the actual constituents of the planetary atmo- 
sphere. Moreover, we also have to relate the real line opacity to 
the effective wide-band absorption assumed here, for which an 
appropriate averaging is required. 

From eq.|7] we can write an expression for the average net 
FIR flux (per unit frequency) over a finite band Av, and define 
an effective opacity through the following: 



AF f 



1 f V2 [ f p,fir(T p ) - ^* ; fir] 



AvJ V| 1 + 3t(v)/4 
[F p ,fir(T p ) - ^+,fir] 
Av 

[F p ,fir(T p ) - F *,fir] 



dv 



(24) 



rV2 



dv 



■ ,-3t(v)/4 
1 



Av 



with Av = Vi — Vj, that is, 
4 



Av 



r- dv 

L 1 + 3t(v 



)/4 



1+3W4' 



- 1 



(25) 



This average is equivalent to the Rosseland mean relevant for 
to the radiation transfer, except that optically thin frequency 
ranges do not cause the flux to diverge. 

Fig.[T2]depicts a specific example, and plots the total effec- 
tive optical depth for two cases with a 20 bands approximation, 
at the high temperature and pressure present on Venus. In the 
first case, only CO2 lines are considered, such that a window 
is left at ~ 70,000A. However, once additional "trace" con- 
stituents are added, the spectral window fills, leaving an opti- 
cally thick spectrum from the optical to the very far IR. For this 
composition, temperature and pressure, we obtain A cut ~ 6000A 
and tfir ~ 10 4 . From fig. [5] we find that the equilibrium tem- 
perature for these values should be T eq ~ 1200K. This is much 
hotter than the real surface temperature on Venus. However, 
the ideal semi-gray model was developed under the assumption 
that the shortwave band is optically thin. This however is not 
the case on Venus because of aerosols and clouds. The finite 
optical depth tsw implies that the solar radiation is deposited at 
different heights, and only a little at the surface itself. An effect 
that lowers the surface temperature. 

In fact, one can think of two limits in which the semi-gray 
model presents a good approximation to reality. The semi-gray 
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Figure 12: The effective band averaged optical depth (as defined in eq. |25( as a 
function of frequency for CO2 at 90 bar and 700K, and for the same atmosphere 
but with added tri-atomic trace elements: 0.014% SO2 and 0.002% H 2 0. The 
line opacity is calculated using the hitran 2004 compilation, and includes ther- 
mal and pressure broadening. Just by itself, the CO2 atmosphere would have 
an optically thin window at ~ 70,()00A, and beyond ~ 250.000A. These would 
allow enough radiation to escape such that the runaway greenhouse would sat- 
urate at relatively low temperatures. With the additional tri-atom molecules 
present on Venus, the windows become optically thick, and A cut moves to even 
shorter wavelengths. 

model corresponds to an optically thin SW band. However, if 
t sw is significant, then SW radiation cannot penetrate the at- 
mosphere and reach the surface before being scattered. This 
implies that the SW energy is deposited at some finite height. 
Up to that height, the atmosphere should be close to isothermal 
while above it, one can assume that the SW band is optically 
thin, implying that the surface temperature should correspond 
to the temperature at this height, and tfir should be the optical 
depth in the FIR of this layer. The case of Venus is therefore 
intermediate. 

8. Summary 

Technically, Simpson's solution assumes implicitly that the 
stellar radiation penetrates the atmosphere, and consequently, 
the greenhouse effect is equivalent to an atmosphere illuminated 
in the IR but from below. The result is that the radiation energy 
density at the bottom of the atmosphere grows linearly with the 
optical depth (in the IR), and the temperature increases like the 
the forth root of t hr . We have shown that this behavior breaks 
down even under the semi-gray approximation. 

Using the full semi-gray model it was shown that two pri- 
mary mechanisms are responsible for the saturation of the run- 
away greenhouse effect. The two cases depend on the value of 
A cut . Unless A cut is small and resides in the optical region, sat- 
uration is achieved by radiating the thermal flux of the planet 
through the SW tail of the thermal distribution. This has an in- 
teresting observational implication: the radiation from such a 
planet should be skewed towards the NIR. 

If A cut is very small, then the saturation of the greenhouse 
could in principle be at a very high temperature, but still less 



than T t . This is far from being realistic with the molecules 
known to exist in planets. The fact that the spectrum is opti- 
cally thin beyond the rotational absorption lines limits the tem- 
perature and can even give rise to an inverse-greenhouse effect. 
Although windows in the FIR for smaller wavelengths could 
cause saturation at a lower temperature, the presence of more 
than one spectroscopically active specie with 3 or more atoms 
can cover the spectrum efficiently. 

Together, the mechanisms imply that a runaway planet like 
Venus cannot have a surface temperature higher than about 1200 
to 1300K. This is significantly lower than some values found in 
the literature. If convection arises, it could reduce the maximal 
temperature even more. 



Regretfully, our long time collaborator and close friend, 
Professor Rainer Wehrse, who actively participated in the last 
draft of this paper, deceased 8.12.2009. 
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Appendix: A simple analytic estimate of the saturation tem- 
perature 



The equation describing the semi-gray model, eq. 12 does 



not yield to an analytical solution. Nevertheless, an approxi- 
mate analytical solution can be obtained when several approxi- 
mations are valid. In particular, if A cul resides in the Rayleigh- 
Jeans long wavelength tail of the distribution of the stellar ra- 
diation, but it also resides in the short wavelength cutoff of the 
thermal emission, then the Planck distribution can be approxi- 
mated and an analytical solution is obtained. This can be used 
to demonstrate the saturation more vividly. 

Specifically, the short wavelength component of the plane- 
tary emissivity can be described by the short wavelength limit 
of the Planck distribution, while the long wavelength compo- 
nent can be described as an integral over the whole Planck dis- 
tribution minus the short wave component. For the incoming 
stellar component, we can do the opposite. The long wave- 
length contribution is described with the long wavelength limit 
of the Planck distribution while the shortwave contribution is 
described by the whole Planck integral minus the long wave- 
length term. Although 6\ = hc/(kT p A cut ) =s 3 - 30 is gener- 
ally large, which allows the above approximate description of 
the planetary component, the appropriate dimensionless num- 
ber for the shortwave stellar component, 62 = hc/(kTi,A cut ) =s 
0.5 - 2, is not always small and the above description can faiQ 

In the short wavelength, we assume that 



2hc I he 



(26) 



2 Note however that the peak of the Planck distribution is at he/ A ~ fewx W, 
therefore a ratio of unity is still on the low side. 



In the long wavelength, we can write 
2 



Ba,fir(T) 



A 1 



kT. 



The required integrals are therefore given by: 



and 
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(29) 



Substituting the approximations of the Planck function for the 



two extreme limits into eq. 12 yields, after some algebra and 
replacing the planet's temperature T p by T sat , that 
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(30) 



where T eq is the equilibrium surface temperature assuming no 
atmosphere. Furthermore, we denote 



a = 



JC I kT eq A cut 



15 



and 
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With these definitions, the equation becomes 



expl - lx(l + 3x). 



(31) 



(32) 



(33) 



We expand the equation around l/a - exp(l/^o), namely, we 
write x - xo+dx, and obtain for the saturation temperature that: 



he 
kA cut 
1 

In a 



(34) 



-i-Mr 5 - 



3 

+ 1 — + 
In a 



a/3 ^ 
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The assumption of a small 6x implies that the first term in the 
square brackets is much larger than the other term (and it can 
be used as a "zeroth order" approximation) . 

A comparison between the numerical and the approximate 
analytic model for the saturation temperature is shown in fig. 



13 The figure reveals that the analytic approximation (calcu- 



lated for the case a — 0) under-predicts the correct T sat by only 
20K to 40K. Most of the discrepancy can be explained through 
the fact that the analytical approximation neglected the stellar 
contribution to the FIR, which although small, does contribute 
some warming. 
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Figure 13: A comparison between the accurate numerical solution of the radia- 
tive transfer equation (eq. |12) and the simple analytical approximation. The pa- 
rameters are those of Venus but with an albedo of a = 0. In addition tfir = 250. 



Appendix B: The full numerical solution 

The main results of the paper are borne from the approxi- 



mate semi-gray treatment, described by eq. 12 However, since 
it is only an adequate description when the SW is optically 
thin, and it ignores thermal redistribution in the atmosphere, 
it is worthwhile to compare the results to the exact solution of 
the semi-gray problem obtained from a full numerical analysis. 
Here we describe the algorithm. 

Basic equations 

In the two stream approximation, the transfer equation for 
the intensities of rays in the upward direction I + (z, A) is 



= - X (z, A)I + (z, A) + X (z, A)S (z, A), 

dz 



(35) 



while the equation for the downward direction 7_(z, A) is given 
by 



- = -xiz,A)I-{z,X) +x(z,A)S(z,A). 
dz 



(36) 



Here the extinction coefficient x(z, A) is the sum of the absorp- 
tion coefficient k(z, A) and the scattering coefficient o~(z, A) at 
height z and wavelength A. T(z) is the temperature at height z. 
Also, the source function S is given by 

S(z,A) = (1 - e(z, A))J(z, A) + e(z, A)B(T(z), A), (37) 

with e(z, A) = k(z, A)/(k(z, A) + <r(z, A)) and B(A, T) being the 
Planck function. 

The mean intensity J is given by 



J(z,A) = (h(z,A) + Uz,A)) 12. 



(38) 



1+ has to fulfill the bottom boundary condition / + (0, A) = 
a/_(0, A) + I + . sm -(A) where a is the albedo of the surface and 
^+,sur(^) is the specific intensity emerging from the surface. The 
boundary condition at the top of the atmosphere is 7_(z max , A) = 
IJA). 



The discretized equations and their solution 

We discretize the monochromatic specific intensities and in- 
troduce the monochromatic intensity vector. 

1(A) = (7 + (zi , A), 7_(z 2) A), 7 + (z 2 , A),..., I + (z N , A)) T (39) 

We describe the radiative transfer by means of transmission 
and reflection coefficients (cf. Shaviv & Wehrse, 1991, Peraiah, 
1984, Wehrse, Baschek & Waldenfels, 2008) 



4 exp(- T ) 



and 



(1 + Ve) 2 " (1 " Vi) 2 exp(-2 Vir) 

(l-6)(l-exp(-2yir)) 
(1+ Vi) 2 -d- Vi) 2 exp(-2Vir)' 



(40) 



(41) 



which are functions of the (wavelength dependent) optical depth 
in extinction t(A) and e(A) - k(A)/(k(A) + cr(A)). 

The transfer equation for the intensity vector then reads 



I = MI + OB + I, 



he 



I = (l-M)- 1 (OB + I 6c ) 



with 
M(A) 



(42) 



(43) 



(44) 
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and a being the albedo of the surface, and 
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a wavelength dependent (2N +l,N) matrix. It has the elements 

ai (A) = (l-t i (A) + r i (Ay)/T i (A)-t i (A) (46) 
bi(A) = (( i W-/- i (i)-l)/T,« + (l-r,W). (47) 

The transfer equation further has 



and 



B = (Z?(zi , A), B(z 2 , A),..., B(zn, A)) 



i bc = (i s «n(A),o,o,...,oy 



(48) 



(49) 
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The vector of the mean intensities is now given by 
J = *FI + l bc /2, 
with the wavelength independent N, (2N +1) matrix 



(50) 



T =2 



1 

110 

110 

1 1 



^\ 















110 
11 
1 1 J 



(51) 



such that the energy equation reads in this nomenclature 

> 

k(A) [B(A) - "F (1 - M^r 1 0(^)B(/l)] dA 
K(A)V(l-M(A)T l l bc (A)dA. 



£ 



(52) 



Boundary conditions 

The lower boundary condition: We assume that the energy 
exchange between surface and the atmosphere is in equilibrium, 
namely: 



oTL + aUO) + h sm (T sw - T air ) = /_(0), 



(53) 



where /(z)_ and /(z)+ are the downward and upward specific 
intensities respectively. /i sul is the heat transfer coefficient be- 
tween the surface and the air through the surface boundary layer. 

The top boundary condition: The top boundary condition is 
the incident stellar radiation, namely 



I— (£ — Zatm) — ^given.stellar- 



(54) 



The outgoing flux is not specified and the converged solution 
should yield an equilibrium, namely vanishing net energy ab- 
sorption by the planetary atmosphere. 
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Figure 14: The deviation from LTE for several optical depths in the particular 
case of A C[U = 13640A. Despite the very large optical depth, the deviations 
from LTE are non-negligible. The leakage of the thermal surface radiation into 
the short-wavelength range varies with height (optical depth), as does the effect 
of the insolation in the long-wavelength range. 



Deviations from LTE 

The numerical solution allows us to include deviations from 



LTE as is demonstrated in fig. 14 in the case of tsw ^ 1 and 
tfir » 1. Plotted is the function iff = ( J - B)/( J + B) where 
J - § MO/47T. In the visible range J » B so that if/ * 1 . If the 
thermal radiation were in equilibrium, then \fj — 0. However, 
we find major deviations close to A ~ A cut . The surprise is that 
even at depths of tfir ~ 200 we notice significant deviations. 

Clearly, if there are windows in the IR (or VIS), then similar 
deviations are expected around the discontinuities in the optical 
depth. 
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